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We explore theoretically how to tune the dynamics of a quasi one-dimensional harmonically 
trapped Bose-Einstein condensate (BEC) due to an additional red- and blue-detuned Hermite- 
Gaussian dimple trap (HGdT). To this end we study a BEC in a highly non-equilibrium state, 
which is not possible in a traditional harmonically conhned trap. Our system is modeled by a time- 
dependent Gross-Pitaevskii equation, which is numerically solved by the Crank-Nicolson method in 
both imaginary and real time. For equilibrium, we obtain a condensate with two bumps/dips which 
are induced by the chosen TEMqi mode for the red/blue-detuned HGdT, respectively. Afterwards, 
in time-of-flight dynamics, we examine the adherence/decay of the two bumps/dips in the conden¬ 
sate, which are induced by the still present red/blue-detuned HGdT, respectively. On the other 
hand, once the red/blue HGdT potential is switched off, shock-waves or bi-trains of gray/dark pair- 
solitons are created. During this process it is found that the generation of gray/dark pair-solitons 
bi-trains are generic phenomena of collisions of moderately/fully fragmented BEC. Additionally, 
it turns out that the special shape of generated solitons In the harmonically trapped BEC hrmly 
depends upon the geometry of the HGdT. 

PACS numbers: 03.75.Lm, 67.85.Hj, 67.85.Jk, 05.45.Yv, 67.85.De, 03.75.-b, 37.10.De 


I. INTRODUCTION 

The development of laser cooling has accelerated a 
tremendous interest in the confinement and manipulation 
of cold atoms. In particular, using optical dipole traps 
generated by red/blue-detuned resonance laser light has 
become a versatile tool for the manipulation of atoms. 
Recently, optical dipole traps created by a red-detuned 
laser beam have become common, experimentally [1- 
12] as well as theoretically [13-19], they are known as 
a “tweezer” or “dimple-trap” provided that the trap is 
quite sharp. Red-detuned dimple traps have become im¬ 
portant tools for BEC production [1, 4, 11], transport of 
a BEC over long distances [5], and formation of shock- 
waves in harmonic plus a dimple trap ]19]. Blue-detuned 
optical dipole traps, instead, are mostly used as a re¬ 
pulsive obstacle for atoms [19-23]. The HG laser beams 
are higher order solutions of the paraxial wave equation 
with rectangular symmetry about their axes of propaga¬ 
tion [24, 25]. Due to their enormous application areas, 
there have been several attempts to develop such higher 
order beam modes [26-28]. For example, the resonator 
of a laser is manipulated such that the beam is emitted 
in a desired beam mode structure [29], or transforms a 
general Gaussian laser-beam with interferometric meth¬ 
ods into the desired modes [30-32] . These interferometric 
methods are typically based on the addition or subtrac¬ 
tion of different scalar laser beam modes [33]. To switch 
between different modes, a more flexible way is to use 
a spatial light modulator to generate the desired higher 
order laser modes [34, 35[. It is already known that HG 
laser modes possess interesting properties, however, to 


the best of our knowledge, so far only two experimen¬ 
tal papers have been published about the confinement of 
atoms in higher order optical dipole traps [36, 37[. 

In this paper, we consider a theoretical analysis of 
a quasi one-dimensional (ID) Bose-Einstein conden¬ 
sate confined by both a harmonic trap and a Hermite- 
Gaussian dimple trap (HGdT). The red/blue-detuned 
HGdT can be generated by using the Hermite-Gaussian 
(HG) laser beam. The mean-field description of the one¬ 
dimensional macroscopic BEG wave function is based 
upon the Gross-Pitaevskii equation [38-40[. A truly 
ID mean-field regime, also known as Tonks-Girardeau 
regime, requires transverse dimensions of the trap on the 
order of or less than the atomic s-wave scattering length 
[41-43[. In contrast, the quasi one-dimensional regime of 
the Gross-Pitaevskii equation holds when the transverse 
dimension of the trap is larger than or of the order of 
the s-wave scattering length and much smaller than the 
longitudinal extension [44-48[ . Here we focus our atten¬ 
tion to a quasi one-dimensional Gross-Pitaevskii equa¬ 
tion (IDGPE). This regime is quite interesting, as it is 
well-known to feature bright solitons for attractive s-wave 
scattering lengths [49-52[, gray/dark solitons for repul¬ 
sive s-wave scattering lengths [53-57[ or the formation of 
shock waves in a BEG [58, 59[. 

With this, we organize our paper as follows. We derive 
the underlying quasi one-dimensional Gross-Pitaevskii 
equation (IDGPE) in Sec. H, where, we also outline the 
system geometry and relate our simulation parameters to 
tunable experimental parameters. Afterwards in Sec. HI, 
for the equilibrium properties of the system, we compare 
a Thomas-Fermi approximate solution with numerical re- 
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suits and show that the HGdT imprint upon the con¬ 
densate wave function strongly depends upon whether 
the HGdT is red or blue-detuned. Later, in Sec. IV we 
assume that the magnetic trap is switched off and we 
determine the time-of-flight (TOF) dynamics of the con¬ 
densate wave function, when the HGdT is still present. 
On the one hand we obtain that for red-detuning the 
HGdT imprint does not decay, but for the blue-detuning 
the HGdT imprint decreases during TOF. On the other 
hand, we discuss in detail how the collision of the con¬ 
densate with the HGdT potential during the non-ballistic 
expansion leads to characteristic matter-wave stripes. 
In Sec. V, we investigate instead matter-wave interfer¬ 
ences in form of the formation of shock-waves/gray(dark) 
pair-soliton bi-trains in the harmonic trap, after hav¬ 
ing switched off the red/blue-detuned HGdT potential. 
There, we also find out that the generation of gray/dark 
pair-solitons bi-trains represents a generic phenomenon 
of collisions of moderately/fully fragmented BEG, which 
strongly depends upon the equilibrium values of the 
red/blue-detuned HGdT depth. Finally, Sec. VI provides 
a summary and conclusions. 


II. MODIFIED QUASI ID MODEL 

We consider a one-component BEG with time- 
dependent two-particle interactions described by the 
three-dimensional GPE 

I! f }^/>(r,t), 

where ^ (r, t) denotes the macroscopic condensate wave 
function for the ®^Rb BEG with the spatial coordinates 
r = (a;, y, z). Here me stands for the mass of the 
®^Rb atom, G® = aB/niB represents the three- 

dimensional ®^Rb coupling constant, where Nb = 20 x 
10^ denotes the number of ®^Rb atoms, and the s-wave 
scattering length is ub = 94.7 ao with the Bohr radius ag. 
Furthermore, V (r) = itibuj^z'^/ 2+mBUj'^ -|- /2 de¬ 

scribes a three-dimensional harmonic confinement, which 
has rotational symmetry with respect to the 2 ;-axis. The 
oscillator lengths for experimental parameters are Iz = 
^hjmB'j^z = 4.12/rm and Zr = i/S/meWr = 0.84^m 
for the trap frequencies = 27r x 6.8 Hz and Wr = 
27r X 160 Hz, respectively. 

An additional three-dimensional narrow Hermite- 
Gaussian laser beam polarizes the neutral atoms 
which yields the HGdT potential I/® = t/oI„m (r). 
Within the rotating-wave approximation its amplitude 
is Go = Sttc^F/ (2wiA) [24, 25, 60, 61], where F = 
|< e|d|g >'^ ui\/ (jiTreohc^) denotes the damping rate due 
to energy loss via radiation, which is detected by the 
dipole matrix element between ground and excited state. 


Furthermore, A = w — wa represents the detuning of the 
laser, here w is the laser frequency and wa stands for 
the atomic frequency. And !„„ (r) describes the inten¬ 
sity profile of the TEM„,^ Hermite-Gaussian laser beam, 
which is assumed to propagate in y-direction and is de¬ 
termined via 

Inm (r) = 2"+’"n!m!7r-^" ((lufl)) 

_[ 2x2 2^2 1 

e Lw 2 (h) w 2 (a)J 

^ WxMIUG) ’ 

with P = f f Inm (r) dxdz being the normalization con¬ 
stant. Furthermore = Wq^/z (l + ZZ^/zzL/z) de¬ 

notes the Gaussian beam radius in the x- and z-direction, 
where the intensity decreases to 1/e^ of its peak value, 
Vrx/z = represents the so-called Rayleigh- 

lengths, which are defined as the distance from the fo¬ 
cus Wqx/z position where the beam radius increases by a 
factor of y/2 [24] . Here Hn (q) and Hm (q) are Hermite 
polynomials of order n and m in x- and z-directions, 
respectively. In the following we restrict ourselves to a 
HGdT potential for a BEG, which is based on a Hermite- 
Gaussian TEMoi laser beam mode and thus carries a 
dark spot in the center of the profile: 


loi (r) = 


8Pz^ 


TrWx (y) (y) 


w2(a) w2(a) 


( 3 ) 


For the TEMqi laser beam, we use the width along the 
cc-axis Wqx = 1.1 yxn and along the z-axis Wqz = 3.2 /rm. 
Therefore, the Rayleigh lengths for the red-detuned laser 
light with A = 840 nm [12] yield y^x = 4.526 p,m and 
2/Rz = 38.29 /im as well as for the blue detuned laser light 
with A = 772 nm [62] we get yBx = 4.92/im and y^z = 
41.6/im. With keeping in mind the fact ?/rx/z 3> l^, we 
can approximate the widths of the HG laser beam in 
X- and z-direction according to W^/z(]j) « Wox/z- This 
simplifies the HGdT potential to 



( 4 ) 


As we have an effective one-dimensional setting due to 
Wz ^ Wr, which implies Iz > h, and j/r,x/z ^ h, we fac¬ 
torize the BEG wave-function via t) = Z)(/(r_L, t) 

with rx = {x, y) and 


(/(rx,t) 


9?^ 


ydrrlx 


( 5 ) 


We follow Ref. [63] and integrate out the two transver¬ 
sal dimensions of the three-dimensional GPE. After some 
algebra, the resulting quasi one-dimensional GPE reads 

-HCb I! tlj{z,t) \\‘^}'tp{z,t), 
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FIG. 1: (Color online) Outer Thomas-Fermi radius Rtfi (red solid), central Thomas-Fermi radius Rtf2 (blue dashed) and 
inner Thomas-Fermi radius Rtf3 (dotted black) as a function of HGdT depth U in dimensionless units for the coupling constant 
value Gb = 11435.9. Below the critical value HGdT depth U < Uc, the BEG is moderately (Mod.) fragmented (Frag.), with one 
exceptional case at 17 = 0, where the BEG is completely confined in one-dimensional harmonic trap. The BEG fully fragments 
into three parts above the critical value Uc ~ 3079, as can be seen in the specific regional inset density plots. 


where V {z) = rfi'Qio\z^ j2 represents an effective one¬ 
dimensional harmonic potential from the MOT, and the 
one-dimensional two-particle interaction strength turns 
out to be 

Gb = . (7) 

Furthermore, the one-dimensional HGdT depth results 
in 


In order to make the IDGPE in (6) dimensionless, we 
introduce the dimensionless time as t = ujzt, the dimen¬ 
sionless coordinate z = zjlz-, and the dimensionless wave 
function ip = With this Eq. (6) can be written in 

dimensionless form 

d ~ f 1 1 T~r~2 

= + r(9) 

+Gb II (5,0 lOj 0 (5,0 , 


here Gb = 2N-QL>j^a-Q/u}zh^ and U = 
SUoPy are the dimension¬ 

less two-particle coupling strength and the dimensionless 
HGdT depth, respectively. The above mentioned exper¬ 
imental values yield the dimensionless Rb-Rb coupling 
constant Gb = 11435.9 and a = Wot./ (V^lz) — 0.548 
represents the ratio of the width of the HGdT potential 
and the harmonic oscillator length along the 2 ;-axis. 
Furthermore, the typical depth of dipole potential traps 
ranges from micro-kelvin to nano-kelvin [64, 65], which 
yields U to be of the order of up to few thousands. From 
here on, we will drop all tildes for simplicity. 


III. STATIONARY CONDENSATE WAVE 
FUNCTION 

In order to determine the equilibrium properties of the 
red/blue-detuned HGdT potential imprint on the con¬ 
densate wave function, we solve the quasi IDGPE (9) 
numerically by using the split-operator method in imag¬ 
inary time [66-69]. The HGdT imprint induces two 
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FIG. 2: (Color online) (a) Height/depth and (b) width of red/blue-detuned HGdT induced bumps/dips according to Eq. (13) 
and W = 2zMa.Ki, respectively, versus red/blue-detuned HGdT depth U in dimensionless units for the experimental BEG coupling 
constant Gb = 11435.9 determined numerically by solving the quasi IDGPE (9) in dimensionless imaginary time (blue circles) 
and analytically (red stars) from the BEG TF wave function (10). 


bumps/dips at the center of the BEG density for 
red/blue-detuned HGdT as shown in the insets of Fig. 1. 
For stronger red-detuned HGdT depth values the two 
bumps increase further, but for stronger blue-detuned 
HGdT depth the two dips in the BEG density get deeper 
and deeper until the BEG fragments into three parts as 
shown in the inset of Fig. 1. To investigate this scenario 
in more detail, we argue that, due to Gb S> 1, the TF ap¬ 
proximation is valid, as the inequality Eint,pot/Ekin 1 
holds within the whole region of interest for the HGdT 
depth U [19]. 

Therefore we perform for the condensate wave function 
the ansatz Tp{z, t) = insert it into the modified 

quasi IDGPE (9), and neglect the kinetic energy term, 
yielding the density profile 





( 10 ) 


In view of the normalization 2 ' II V’W f dz = 1, 

which fixes the chemical potential /r, we determine the 
Thomas-Fermi radii Rtf from the condition that the con¬ 
densate wave function vanishes: 


/,= ^+C/R^^e-’^. (11) 

As can be read off from the inset of the Fig. 1 the 
number of solutions of Eq. (10) changes for increasing 
red/blue-detuned HGdT depth U. In the case, when U 
is smaller than Uc, Eq. (10) defines only the BEG cloud 
radius Rtfi- But for the case U > Uc, the blue-detuned 
HGdT drills two holes at the center of the ®^Rb conden¬ 
sate, so the BEG fragments into three parts as shown 


in Fig. 1. Thus, we have then, apart from the outer 
condensate radius Rtfij also two inner condensate radii 
Rtf 2 and Rtfs- With this the normalization condition 
2 fo II V' (z) \\^ dz +2 II V' (z) f dz = l yields 


M (R-tfi — Rtf2 + Rtfs) ~ g (Btfi ~ B-tf2 + Btfs) (^2) 


R' 


TFie 


— RTF2e 


- Rtfsg 




EU(^]-EU(^]+ Erf ( 


V ^ 


- 


V 




\ ^ 


- 


where Erf(y) = e ^ dx denotes the error function. 

In case oi U < Uc, the BEG’s inner two radii Rtf 2 and 
Rtfs vanish and the BEG outer radius is approximated 
via Rtfi ~ due to Eq. (11). Thus, for U < Uc 
the BEG chemical potential is determined explicitly from 
(12): ^ « 3^/^/2’^/^ (2 Gb + . Provided that 

U > Uc, two inner cloud radii Rtf 2 and Rtfs have to be 
taken into account according to Fig. 1. We observe that 
the Thomas-Fermi value of the critical red/blue-detuned 
HGdT depth Uc « 3079 is close to the numerical one 
Uc ~ 3090. Figure 1 also shows the resulting outer and 
inner Thomas-Fermi radius as a function of the red/blue- 
detuned HGdT depth U. Here, the two inner radii behave 
symmetric, e.g., for U > 3079 the Rtf 2 is increasing and 
Rtfs is decreasing correspondingly, however after U > 
4500, they both become approximately constant as shown 
in Fig. 1. We also read off that Rtfi ~ remains 
approximately constant for U > Uc, so we conclude that 
the chemical potential /r is then sealed to its critical value 
He « 341.28. 

In the perspective of a quantitative comparison be¬ 
tween the analytical and the numerical calculation, we 
characterize the red/blue-detuned HGdT induced im¬ 
print upon the condensate wave function by the follow- 
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FIG. 3: (Color online) Height/depth versus dimensionless time during TOF for different (a) red-detuned and (b) blue-detuned 
HGdT depths U in dimensionless units. 


ing two quantities. The first one is the red/blue-detuned IV. TIME-OF-FLIGHT DYNAMICS OF 

HGdT induced imprint height/depth RED/BLUE-DETUNED HGDT INDUCED 

IMPRINT 


r Max(|| ||2) - 11 ^-(0) ||2 U<0 

HD = I (13) 

i II vI/(0) f-Min(|| vI/(^) ||2^„) U>0 


and the second one is the red/blue-detuned HGdT in¬ 
duced imprint width W = 2zMax, where ZMax denotes the 
coordinate of maximal density. To find out a one-to-one 
resemblance between analytical and numerical calcula¬ 
tion of HD and W, we determine the solution of the di¬ 
mensionless IDGPE (9) and compare it with the TF solu¬ 
tion of Eq. (10), as shown in Fig. 2. The case U = 0, i.e., 
when the HGdT potential is switched off, corresponds to 
a BEG in a quasi one-dimensional harmonic trap. Fur¬ 
thermore, in the range U < f7c we observe that the 
red/blue-detuned HGdT induced imprint height/depth 
changes linearly with the optical dipole trap depth U ac¬ 
cording to 


HD 




U 7^ 0 , a 7^ 0 , 


(14) 


where T = xe^ abbreviates the productlog function. In 
case of [/ > Uc the blue-detuned HGdT induced im¬ 
print depth yields the constant value HDc « 0.0296 as 
follows from Eq. (14), which slightly deviates from the 
corresponding numerical value HDc = 0.027. Similarly, 
the red/blue-detuned HGdT induced imprint width fol¬ 
lows from W = ^OL\/\ — T (—e/2|[/|) according to the 
TF approximation, which reduces at the critical blue- 
detuned optical dipole depth to Wc « 3.60 whereas the 
corresponding numerical value is Wc ~ 3.64, as shown in 
Fig. 2. 


The time-of-flight (TOF) expansion has been used to 
measure various BEG properties since the field’s incep¬ 
tion. By suddenly turning off the magnetic trap, when 
the HGdT is still present, the atom cloud is allowed to ex¬ 
pand in all directions. This expansion proceeds according 
to the momenta of the atoms at the initial time t = 0 and 
an additional tiny acceleration results from inter-particle 
interactions. The red-detuned HGdT induced two bumps 
remain approximately constant during the temporal evo¬ 
lution as shown in Fig. 3(a). But the blue-detuned HGdT 
induced two dips at the center of the condensate start 
decaying with a characteristic time scale after having 
switched off the trap as shown in Fig. 3(b). Further¬ 
more, the dips of the HGdT induced imprint start de¬ 
caying faster with increasing blue-detuned HGdT depth 
for smaller time as shown in Fig. 3(b). Note that the rel¬ 
ative speed of the bumps or dips from each other turns 
out to vanish. 

Furthermore, we investigate in detail the possible oc¬ 
currence of matter-wave stripes at the top of the con¬ 
densate during the non-ballistic expansion of the moder¬ 
ately/fully fragmented BEG cloud by plotting the den¬ 
sity distribution of the released cloud correspondingly 
as shown in Fig. 4. According to Fig. 4(a) we do 
not observe any particular structure for the red-detuned 
HGdT induced imprint, but Fig. 4(b-d) shows for the 
blue-detuned HGdT induced imprint that characteristic 
matter-wave stripes occur, which are generated, while the 
freely expanding BEG collides with the HGdT potential. 
For small blue-detuned HGdT depth, the generation of 
matter-wave stripes can be seen at later time, as com- 
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FIG. 4: (Color online) Time-of-flight evolution of de- 
pleted density || ipiz,t) ||dd = II ||y^o - 

II Ip {z,t} Wij^Q from solving the modified quasi IDGPE equa¬ 
tion for different values of |17|: (a) U = —500, (b) U = 1500, 
(c) U = 2500, and (d) U = 3500 in dimensionless units. 


pared to higher blue-detuned HGdT depth, as shown in 
Fig. 4(b-d). The matter-wave stripes are directly vis¬ 
ible for U < Uc, as can be explained as follows. In 
Fig. 4(b), the height of the blue-detuned HGdT induced 
two dips is smaller as compared to Fig. 4(c), therefore 
they need more time to drill a hole in the condensate dur¬ 
ing TOF. For the HGdT potential depth U = 1500, the 
BEG fragments into three parts at the dimensionless time 
t = 2.4, afterwards the three fragmented condensates 
start to interact as separate identities with the HGdT 
potential, which leads to the formation of characteristic 
matter-wave stripes. The similar phenomenon happens 
in Fig. 4(c), but in this example the initial HGdT poten¬ 
tial depth U = 2500 is larger than the previous one in 
Fig. 4(b), so the BEG becomes fragmented at the earlier 
time t = 1.3. In the example of Fig. 4(d), when U > Uc, 
the BEG is already initially, i.e. at time t = 0, frag¬ 
mented according to Fig. 1. Therefore the matter-wave 
stripes can be seen just after t > 0, but the stripes are 
not as visible as in the two previous cases. 


V. SHOCK-WAVES AND GRAY/DARK 
PAIR-SOLITONS BI-TRAINS 

In this section, we show that matter-wave self¬ 
interferences emerge, once the red/blue-detuned HGdT 
potential is suddenly switched off, within the remaining 
harmonic confinement, as this leads to shock waves and 
gray/dark pair-solitons bi-trains, respectively as shown 
in Fig. 5(a-d). A shock-wave is a special kind of propa¬ 
gating disturbance in the BEG, whose amplitude, unlike 
for solitons, decreases relatively quickly with large dis¬ 
tance. Furthermore, gray/dark solitons have a charac¬ 
teristic property that they can pass through one another 
without any change of shape, amplitude, or speed. We 
can see from Fig. 5(b-d) that the pair-solitons bi-trains 
do, indeed, pass through one another and that they are 
reflected from the end of the trapping potential. 

Once the red/blue-detuned HGdT potential is 
switched off, the system quasi-instantaneously adjusts 
its energy to the new equilibrium, paving the way for 
the creation of shock-waves and bi-trains of gray/dark 
pair-solitons, respectively. The total normalized energy 
E(t) = E(t)/Max(E(t)) as shown in Fig. 5(i-l), changes 
quite quickly from its initial value to a new equilib¬ 
rium value, thus generating the shock-waves or the pair- 
solitons bi-trains. For an initial red- and blue-detuned 
HGdT depth, we observe that two excitations of the 
condensate are created at the position of the red/blue- 
detuned HGdT potential, which travel in the opposite di¬ 
rection with the same center-of-mass speed, are reflected 
back from the harmonic trap boundaries, and then col¬ 
lide at the red/blue-detuned HGdT potential position as 
shown in Fig. 5(a-d). 

We have performed calculations for different red- 
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FIG. 5: (Color online) Coherent matter-waves evolution in the BEC density, when red/blue-detuned HGdT potential is switched 
off, versus time and position for different values of (a) U = —500, (b) U = 1500, (c) U = 2500, and (d) U = 3500. Graphs 
(e-h) show the BEC density at time t = 5.3, and in graphs (i-1) we plotted the corresponding energies E{t) = E{t)/Max.{E{t)) 
versus time for the HGdT depths U of the graphs (a-d) in dimensionless units. 


detuned HGdT potential depths U < 0 and in all cases 
we observe the formation of the shock-wave structures 
as shown in Fig. 5(a). The density of atoms around the 
shocks is mostly enhanced in comparison with the den¬ 
sity far away from these perturbations. And for the blue- 
detuned HGdT potential trap, we detect gray/dark pair- 
solitons bi-trains, traveling in opposite directions with 
the same speed as shown in Fig. 5(b-d). The creation 
of these calibrated gray/dark pair-solitonic bi-trains are 
generic collision phenomena of moderately/fully frag¬ 
mented BEG, which is strongly depending upon the equi¬ 
librium values of the red/blue-detuned HGdT potential 
depth, respectively. 

The dynamics of one gray/dark soliton in a BEG cloud 
is well described by z = —(l/2)9I4xt/9z, where V/xt is 
the dimensionless confining potential and z denotes the 
position of the gray/dark soliton. In the case of har¬ 
monic confinement with a potential V/xt = 2^/2 the so¬ 
lution of this evolution equation leads to an oscillation 
of the soliton described by z{t) — Rtfi sin (f/-\/2). Thus 
the frequency of the oscillating soliton and the frequency 
of the dipole oscillation of the Bose-Einstein condensate 
in the trap differ by the factor [70]. In our system, 
pair-solitons bi-trains generally oscillate with the average 
frequency fl = 27r x 4.80 Hz irrespective of the sign and 
the size of U as shown in Fig. 5. With this, we get the 


ratio n/wz « 0.705, which is quite close to the dimen¬ 
sionless soliton frequency l/-\/2 « 0.707 in a harmonic 
trap as predicted in Ref. [70]. Note that previously the 
generation of solitons was studied theoretically by inves¬ 
tigating the collision of two condensates [71] and experi¬ 
mentally for different quasi one-dimensional trap geome¬ 
tries [72, 73]. Although in the latter experiments only 
one potential maximum occurs instead of two as in our 
work, so there single solitons and here pairs of solitons 
are observed, the basic physics is the same. 

We also observe an intriguing substructure of each soli¬ 
ton, which we call pair-soliton. Normally, we find that 
there are solitons which always move in pairs, and the 
mean distance between each other is less than the neigh¬ 
boring solitons as shown in Fig. 5(b-d) and Fig. 5(f-h). 
Numerically, we have observed that the averaged distance 
between pair-solitons is less for dark solitons as com¬ 
pared to the gray solitons as shown in Fig. 5(f-h). We 
also observe that, in general, a minimal time of about 
4.6ms is required to generate shock-waves/pair-solitons 
bi-trains as shown in Fig. 5. The number of shock-waves 
is not effected by the red/blue-detuned HGdT poten¬ 
tial depth, but the number of interference fringes in¬ 
creases. On the other hand, we observe that the number 
of gray/dark pair-solitons depends on the depth of the 
red/blue-detuned HGdT potential as shown in Fig. 6, 












































FIG. 6: (Color online) Number of pair-solitons Nsoi in each 
solitonic train versus HGdT potential depth U in dimension¬ 
less units. 

the highest number of pair-solitons in every train is 7. 
For the blue-detuned HGdT depth U <Uc, the number 
of pair-solitons grows linearly in the condensate and after 
the critical value U^, the number of pair-solitons remains 
approximately constant. 

Note that in case of the collision of two condensates 
in Ref. [71], it turned out that the number of observable 
solitons depends sensitively on the initial phase difference 
of both condensates. Thus, if the two condensates have 
an initial phase difference of 0(7r), the number of solitons 
is even(odd). In our case, we have a single BEG fragment¬ 
ing into three parts, which have the same phase, therefore 
we observe an even number of pair-solitons in the conden¬ 
sate in agreement with Ref. [71]. Indeed, Fig. 6 shows 
the number of pair-solitons in each solitonic train, so the 
total number of pair-solitons in the whole condensate is 
twice as large. But, in Fig. 6 it turns out that the num¬ 
ber of pair-solitons depends crucially on the depth of the 
red/blue-detuned HGdT potential as shown in Fig. 6. 


VI. SUMMARY AND CONCLUSION 

We have developed a simple quasi ID model both an¬ 
alytically and numerically to calculate the statics and 
dynamics of the red/blue-detuned HGdT imprint upon 
the ®^Rb condensate. First of all, we showed a quan¬ 
titative comparison between the Thomas-Fermi approx¬ 
imation and numerical solutions for the underlying ID 
Gross-Pitaevskii equation for the equilibrium properties 
of the proposed system. Later we discussed that the 
HGdT potential imprint upon the condensate wave func¬ 
tion strongly depends upon whether the effective HGdT 
is red- or blue-detuned. With this we found out that the 
HGdT imprint generates two bumps/dips at the center of 
the BEG density of the red/blue-detuned HGdT. Later 
we discussed that the red-detuned HGdT induced two 


bumps did not decay when we switched off the harmonic 
trap but the blue-detuned HGdT induced two dips de¬ 
cay. During the time of flight, we saw the emergence of 
matter-wave stripes at the top of the condensate, which 
arise as the BEG decomposes into a fraction at rest in 
the center and two moving condensates at the borders. 
We have used the quasi one-dimensional time-dependent 
GPE to analyze the creation of gray/dark pair-solitons 
bi-trains within the moderately/fully fragmented BEG, 
which is strongly depending upon the HGdT potential 
depth. The Hermite-Gaussian dimple trap geometry 
maybe more applicable to soliton interferometry rather 
than the Gaussian barrier adopted in Refs. [19, 74], 
because one can shape solitons. Additionally, we also 
showed that the number of pair-solitons in the system is 
depending on the initial HGdT potential depth U. Dur¬ 
ing the generation of pair-solitons it was astonishing to 
And that the special shape of the newly generated soli¬ 
tons in the harmonically trapped BEG is sculptured by 
the external potential and the generation of gray/dark 
pair-solitons bi-trains is a generic phenomenon of colli¬ 
sions of moderately/fully fragmented BEG. With this we 
conclude that it maybe possible in the future to frame 
complex shapes of solitons in the harmonically trapped 
BEG by imposing a unique geometrical configuration for 
the external potential. 

The ability of sculpting a quasi one-dimensional har¬ 
monic trapped Bose-Einstein condensate by a HGdT has 
many exciting prospects. For instance, it can be used to 
generate a truly continuous atom laser, which has many 
applications in atom interferometry [75, 76] . To construct 
such an atom laser one needs a device that continuously 
converts a source of condensed atoms into a laser-like 
beam. In Sec. IV, we saw in the time-of-flight picture 
for the case U > Uc that a BEG reservoir occurs at the 
center of the trap. By suitably tuning the HGdT depth a 
fraction of this fragmented condensate could be coupled 
out, serving as a source for an atomic beam. 
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